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Abstract: We discuss the problem of incorporating recoil effects into the prob- 
abilistic QCD evolution scheme based on the picture of colour dipoles as done in 
recent Monte Carlo programs. Such a scheme correctly describes subleading soft 
contributions to multiplicity distributions. However we find that a simple receipt 
for incorporating recoil effects into the dipole multiplication picture based on the 
ordering of gluon energies conflicts the coUinear factorization and does not lead to 
the correct DGLAP equation. 
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Preface 

We share opinion of our late friend Bo Anderson who once bitterly complained about 
two aspects of scientific publishing that became standard in recent decades: 1) shying 
away from expressing emotions and 2) not making public an outcome of a study that 
produced a negative, undesired results. The present paper is moderately emotional; 
at the same time it is perfectly in line with the second Bo's demand: it reports the 
finding that the authors did not plan, and would rather not want, to arrive at. 

The answer to the question that is posed in this paper suggests that a formulation 
of Monte Carlo event generation based on the multiplication picture of energy ordered 
colour dipoles is plugged by a serious problem. 



- 1 - 



1. Introduction 

Generation of events using Monte Carlo methods is an indispensable tool for plan- 
ning, running and analysing the results of modern high energy experiments. A 
possibility to generate multi-particle production as a Markov chain of successive 
independent parton splittings is based on the general property of factorization of 
collinear singularities which is characteristic for logarithmic field theories (quantum 
field theories with dimensionless coupling). 

Event generation is necessarily an approximate procedure. Strictly speaking, in 
order to predict the cross section dan for production of n partons in a hard interaction 
with the hardness scale Q^, one needs to plug in the QCD matrix element of the (at 
least) n*^ order of the perturbation theory. Instead, one repeats n times the basic 
1 — > 2 sphtting process of the first order in Ug. By so doing one correctly samples 
the major part of the n-parton phase space, though not all of it. It is this part that 
gives the dominant, maximally collinear enhanced contribution to, e.g., the inclusive 
particle distribution, D^^^ = (9(a"log"Q^). This approximation known as leading 
logarithmic (LLA) can be systematically improved by including higher order parton 
splitting processes. Thus, the next-to-leading logarithmic accuracy corresponding 
to the scries a"^^ log" is achieved by incorporating 1 — > 3 parton splittings, etc. 
For inclusive parton distributions this leads to the DGLAP evolution equations [1] 
whose generalisation to multi-particle distributions can be achieved in the spirit of 
jet-calculus [2], with additional account of soft gluon coherence. 

Beyond the one loop "evolution Hamiltonian" it starts to matter how one or- 
ganises the parton chain, that is, what one chooses as "parton evolution time". In 
particular, due to the presence of soft gluons one encounters corrections of the type 
(ttslog^a;)'^ which arc formally subleading (non-collinear) but become explosively 
large in the a; ^ 1 limit (here x is the gluon energy fraction). This eventuality 
can be cured by an appropriate tuning of the evolution time parameter. Namely, 
for time-like parton multiplication (jets) one has to choose the angular ordering [3]. 
This choice corrects the "naive" fluctuation time ordering (dictated by an examina- 
tion of Feynman denominators in a cascading process) by taking into clever account 
collective effects leading to destructive interference in soft gluon radiation (for more 
details see, e.g., [4-8]). 

An alternative way of dealing with soft gluon interference effects for time-like 
parton multiplication (jets) is provided by the "dipole scheme" [4, 9] in which an 
independently radiating parton, 1 ^ 2, is replaced by a colourless dipole formed 
by two partons neighbouring in the colour space, 2 — 3. Gluon radiation off a 
dipole is automatically suppressed at angles exceeding the dipole opening angle thus 
reproducing the angular ordering. Colour suppressed contributions C(l/A^) that 
lie beyond the scope of the dipole approach are numerically small and difficult to 
trigger experimentally [10]. 



- 2 - 



Within the standard coherent parton cascade picture [5,8] it is the small-angle 
multiplication processes populating jets that enjoy full all-order treatment (get "ex- 
ponentiated" ) . The dipole formulation offers a possibility to improve this treatment 
by taking into consideration logarithmically enhanced effects due to multiple emis- 
sion of soft gluons at large angles with respect to jets. Non-coUinear soft gluons 
dominate inter-jet particles flows in various hard processes. They also complicate 
the analysis of the so-called nonglobal QCD observables [11-13], i.e. in observables 
in which recorded radiation is confined in geometrically definite phase space regions. 
It is then interesting to involve these corrections into a Monte Carlo, code based on 
dipole emission, including soft radiation awy from jets. 

In this paper we discuss a dipole based scheme (in the large- A^c approximation) 
which is well suited for deriving improved analytic predictions for observables like 
mean multiplicities, inclusive soft particle spectra and energy flows, correlations, etc., 
that incorporate large-angle soft gluon radiation effects as in [11-13]. At the same 
time, we observe that once one aims at beyond the no-recoil (soft) approximation, 
treating colour dipoles as independently evolving entities is likely to conflict the 
coUinear factorization. It does not lead to the correct DGLAP equation. 

The paper is organised as follows. We start by constructing in Section 2 an 
evolution equation for the generating functional that describes soft gluon content 
of e'^e" annihilation events. No recoil is included. This equation is based on the 
factorization of the multi-soft gluon distribution, in the planar limit, and uses the 
centre-of-mass energy as an "evolution time". Contrary to what is done normally, 
no angular variable enters the evolution time here. The evolution equation in the 
integral form we interpret as a Markov process of successive dipole branchings, which 
interpretation leads to a Monte Carlo chain process for probabilistic event generation. 
As a flrst example, we solve the evolution equation for mean gluon multiplicity 
to analyse subleading corrections due to large-angle soft gluon radiation. Here we 
work in the soft approximation meaning: 1) small energies of emitted gluons and 2) 
neglecting the energy-momentum loss (recoil) by the primaries. 

These two aspects of the "soft gluon approximation" do not necessarily go to- 
gether. In Section 3 we consider an inclusive energy distribution of the leading quark 
in the quasi-elastic limit, {1 — x) <^ 1, where emitted gluons ought to be soft but 
the quark recoil is essential. At the beginning, in the construction of the evolution 
equation no recoil effects were included. They are incorporated later by using the 
recoil strategy that seems natural for a dipole multiplication scheme, see [14]. We 
flnd, however, that the quark fragmentation function so obtained deviates essentially 
from the coUinear resummation result given by the DGLAP evolution equation. This 
is in contrast to what is happening in Monte Carlo schemes based on angular ordered 
time-like cascades where coUinear singularities are correctly resummed [5,15]. 

Discussion of the problem encountered is relegated to the Conclusions section. 
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2. Multiple soft gluons and Monte Carlo 



The method of generating functionals makes it straightforward to generate exclusive 
events using Monte Carlo techniques. We start from the known n-gluon distribution 
in the soft limit and construct the evolution equation for the corresponding generating 
functional. 

The generating functional for a quark-antiquark pair qq plus an ensemble of n 
secondary partons, 7* PaPb gi, ^2, • • • , is 

G{Q,u) = Y, / d^^:^{q,,q2,...,qn)uaU,\[u{qi)\M^\q,,...,qn)\ . (2.1) 

n=0 i=l 

Here pa-, Pb are the final on- mass-shell momenta of the quark and antiquark produced 
by the virtual photon 7* and qi the momenta of the secondaries. The matrix element 
describes production of n + 2 real final state partons and contains any number 
of virtual ones, d^^"'^ is the full phase space factor of the ensemble of n + 2 massless 
particles: 

n n 

g2, . . . , Qn) = {dpa){dpb) l[{dq^) (27r)^5(^)(p, + + J] - Q), (2.2) 

i=l i=l 

where 

{dq) = ^ = go = |q|. 

The "source functions" u are attached to each parton in (|2.1j ). They help to extract 
an arbitrary final state observable once the functional is known. The fully inclusive 
measurement, that is when one allows for production of any number of particles 
with arbitrary momenta, corresponds to setting all u=l. This gives G{Q,u = l) = 1 
corresponding to normalization to the total cross section. In this section we work in 
the no-recoil approximation. 

2.1 Generating functional in the soft limit 

In the large- Ai'c ( "planar" ) approximation the squared matrix element can be approx- 
imated at tree level as 



2 (A-rr'^ 



M'J{qu...,q„) ~ ^ --^\Wab{qi,,...,qi^, 

nl ^ — ' TT 

perm 

(2.3a) 

where Wab is the "multiple antenna" function [4] describing the production proba- 
bility of n soft real gluons g^: 

Wab{qi,q2,---,qn) = 7 w ^^"^''^ . r- (2.3b) 

{Paqi){qiq2) - ■ ■{qnPb) 



2 _ Ar,«, 



a. 
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This result holds in the soft gluon approximation qi <C Pa,Pb- In invariant terms, 

(QiPa), (qiPb) < (PaPb), (2.4) 

or, equivalently, {qi,Pa + Pb) < '2{paPb) = Q"^ giving 

< ~ Efe ~ g/2, (2.5) 

with Ea,Eb and Ui the energies in the centre-of-mass of the "dipole" {ab). In this 
approximation the phase space can be simplified as 

n n 

rf$i:Hgi,...,g„) ^ {dpa){dp,){27ly6\pa + Pb-Q)■l[(d1^) = ^$1? " Hl^?.) • (2-6) 



i=l 



The expression ( |2.3| ) was derived implying strong ordering of gluon energies: 

The answer, however, is symmetric with respect to momenta of participating gluons. 
Therefore, we may ignore the ordering between gluons and set a common upper 
bound E = \Q for gluon energies. 

Extracting the quarks whose momenta in the soft gluon approximation are not 
affected by radiation, we may write 

G{u- Q)c^J uipaHp,) \M^^\' ■ g{pa,Pb; IQ). (2.7) 

The functional g{pa,Pb', E) describes production of any number of gluons with ener- 
gies limited from above by some value E, off the dipole formed by the quarks with 
momenta pa, Pb- By setting E = we obtain the multi-gluon generating functional 
describing e~^e~ annihilation process with s = Q^. 

The generating functional contains proper virtual factors for each n-gluon con- 
tribution. Omitting for the moment the virtual corrections, we have 



/n 
Wabiqi, ■..,qn) Yl 
■n i=l 



ujiduj— as ■ u{qi) ■ ^{E-uJi) 



(2. 



Due to the sjTiimetric structure of the antenna function ( |2.3| ) one may choose Wab{qi ■ ■ ■ qn) 
to represent any energy ordered ensemble of gluons and drop the symmetry factor 
1/n!. 

2.2 Evolution equation 

To obtain an evolution equation for the generating functional we need to exploit the 
structure of the multi-gluon antenna distribution ( p.3| ) and to deduce a recurrence 
relation. To construct a recurrence relation one should select an "evolution variable" 
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to order emitted gluons. In the logic of the collinear approximation this is done by 
ordering the angles of successive gluon emissions, which ordering takes full care of 
the destructive interference contributions in the soft region [3-5,8] and preserves the 
probabilistic parton multiplication picture. In particular, this procedure was applied 
to the expression ( p.8| ) to construct the HERWIG event generator [5,15]. 

However, as was already stated, the distribution ( |2.3| ) is valid in the soft limit 
for arbitrary angles, and one is tempted to use this property in order to lift off the 
collinear approximation. The derivation of the evolution equation avoiding the small 
angle approximation was done in [11-13] where the upper bound E on gluon energies 
in the event centre-of-mass was treated as the "evolution time" parameter. Let us 
recall the corresponding construction employing the energy as an evolution variable. 

2.2.1 Energy as an evolution variable 

We take the ordering in the energy of emitted gluons in the center-of-mass of the 
event (|2.5| ) and deduce an evolution equation for the soft generating functional. The 
evolution equation follows from the exact recurrence relation — the factorization 
property of the multi-gluon antenna function ( p.3b| ) — which reads 



Wabiqi, g2, • • • , qn) = Wabiqi) " Waiiqi, ■ ■ ■ qe-i) Wib{qi+i, ...,qn), (2.9) 

where i is any gluon. To construct the evolution equation one takes i to be the most 
energetic among the gluons. 

Introducing the scaled antenna function 

Wab{q) = uj'^ ■ Wab{q) = 7%- > 6fc = 1 - cos Qik , 

Sag Sfe<7 



which depends only on the angles between partons, we obtain 



uiq)GiPa,q;uj)g{q,pb;uj)-giPa,Pb,uj) 



. (2.10) 



Here the subtraction term takes care of virtual corrections, so that taking all u{q) = 1 
one derive the desired normalization 

g{pa,Pb;E)l^, = l. (2.11) 

No recoil was included up to now; it will be considered later. 

2.2.2 Integral evolution equation and Monte Carlo 

By exponentiating the total one-gluon emission probability we obtain the Sudakov 
form factor which describes the probability that the dipole does not radiate gluons 
with energies up to a given value E: 

S{pa,Pb;E)=expL r^^a.(gf)W?.,(g)%f -Qo)|. (2.12) 
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Here qf' is the invariant transverse momentum of the gluon with respect to the pair 
of quarks: 

and an arbitrary parameter Qo has been introduced as a colhnear cutoff. By con- 
structing the logarithmic derivative of the Sudakov form factor over the maximal 
gluon energy, 



(2.14) 



and plugging this expression into ( |2.1(]| ) we can trade the virtual subtraction term 



for the Sudakov factor. By so doing we arrive at the equivalent equation for the soft 
generating functional in the integral form: 

du f dQ S{pa,Pb;E) 



duo [ dVl 

y[Pa,Pb]E) = S{pa,Pb;E)+ / — / u{q} 



S{pa,Pb;uj) (2.15) 
• asiqf) WM - Qo) ■ G{Pa, q; oo) G{q,p,; uo). 

Iteration of this equation can be interpreted as parton branching which can be re- 
alised as a Monte Carlo (Markov) process. To see this, observe that the iteration 
kernel 

ciP(Pa,P., g; i?) = S{p^^p,,E) ^ _ 

u) An S{pa,Pb;uj) 

can be written as the product of two probability distributions (here the boundary 
gf^ > Qo is imphcit): 



qJ ' 



S{pa,Pb\^) " 47r iaqiqb\J An^^aqUb 



(2.17) 



In the Monte Carlo process one tries to generate emission of a gluon with momentum 
g off a given dipole. The first distribution, dV^^\ provides the gluon energy cj, if the 
transverse momentum bound is satisfied, and the second one, dV^'^\ its direction VLq. 
If the boundary is not satisfied, the dipole does not emit. One has 

dV{pa,Pb,q,E) = l-S{pa,Pb;E), (2.18) 

which shows that the Sudakov factor gives the probability of not emitting a soft 
gluon within the resolution g"^ > Qq. 
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2.3 Mean multiplicity 

By construction the generating functional ( p.l5|) embeds only soft gluon radiation. 
Incorporating double logarithmic (simultaneously soft- and collinear-enhanced) con- 
tributions, it also contains single logarithmic effects due to emission of soft gluons at 
large angles. 

Here we demonstrate an application of the dipole-based evolution to the calcu- 
lation of subleading large-angle soft gluon corrections. In the collinear approach, 
such corrections are treated as due to "multi-jet" configurations (contributing to the 
"coefficient function") rather than jet evolution ("anomalous dimension"), see [8]. 

Consider mean gluon multiplicity. In order to obtain an equation for the multi- 
plicity of secondary partons, one applies to ( |2.1(]| ) the variational derivative over the 
probing function u{k) and integrates over k, while setting u = 1 for all remaining 
probing functions (one-particle inclusive measurement): 



Introducing A/" = 1 + we have [16] 



1 + Ni^ga, E) + iV(e,fe, E) - N{U, E) 



ui=E 

;2.i9) 



^as{kf)WM 



ab) 



(2.20) 



In order to illustrate the difference with the standard approach based on the angu- 
lar ordering we will analyse ( p.20[ ) in the double logarithmic approximation which 
neglects hard parton splittings and recoil effects. For the time being we will ignore 
the running of the coupling; the corresponding corrections will be addressed later. 
Representing the integrand on the r.h.s. of ( p.2(]| ) as 



66 

and using the 1 



1 

6 



Ar(ei) - Ar(0 



+- 



2 symmetry, we obtain 



m) 



66 



m) 



e-6-6 
66 



(2.21) 



We can perform the integration over the azimuth of q around the direction of p^. 



dVL 
Ait 



rf6 



271 



Using the relation 



2vr 6 - 6 
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we get 



^-^1-6 



As a result, 



66 



1 f^-^i 



le-eii 



6 vie -61 



-mi) -mo 



U 6 



(2.22) 



(2.23) 



with the E dependence imphcit. The first contribution we spht into integrals over 
small and large angles, < < ^ and ^ < 6 ^ 2. Introducing the integration 
variable rj = ^i/^ < 1 in the first region and t] = ^/^i in the second, |^ < < 1, the 
equation ( p.23| ) takes the form [16] 



drj as 



v 



lo 1-^ 

This result we represent as a sum of two terms 

drj a 



+ Lp^Mv-'^)-m)]- (2.24) 



EdEATiO 



^ 



(2.25a) 



where 



m= i'?^mvo-m)] + f! ^W{v-'^)-m)]- (2.25b) 



1-^ 



The first integral term on the r.h.s. of ( p.25a|) generates the standard DL anoma- 
lous dimension for the mean multiplicity; the second one constitutes a subleading 
correction. 

2.3.1 The main term 

We have to solve the equation 



V 



EdEM^'^{i;E) 



(2.26) 



with the initial condition 

m;E)\E^^^Q, = 1. 

The solution is a function of a single variable = Ey/2^ (the maximal transverse 
momentum of partons radiated by the dipole). 

The argument of the running coupling in ( p.26|) can be approximated as 
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The correction to this approximation turns out to be neghgible since the logarithmic 
factor In (^2/0 vanishes in the colhnear hmit ^1 — 0: 



oc ||i(/3o«^)-ln| = 0{6?^. 



Thus we obtain the double differential equation 

(dl^) = 2a.(Q5)ArDL(Q^), Q^ = Ev^, (2.27) 

with the initial conditions M{Qq) = 1 and M'{Qq) = 0. Its solution for the fixed 
coupling case reads 

^-W,)^co*(„„|)4(|)\ (2.28) 

where the anomalous dimension 7 is given the correct expression 

- f^dT] I _ 2 /_ 

1 = as r]2^=as-- =^ -f = V2as. 

Jo V 7 

2.3.2 The correction term 

Let us analyse the correction term ( |2.25b|) keeping the contributions of the order of 
a^^'^ while neglecting contributions O (a^). The answer depends on the value of the 
opening angle ^. 

Small opening angles '€i 1). It is the region of small opening angles that gives 
a dominant contribution to the integral determining the DL anomalous dimension 
in (|2.26| ). In this entire region the correction ( p.25b|) does not contribute at the a^^"^ 
level. Indeed, substituting the solution ( p.28|) into the expression ( p.25b|) and setting 
^ ^ in the lower limit of the second integral we get 



= a 



I ^[r,-/^ + r^/^-2] 
2^(l)-^(l + i7)-^(l-|7) 



(2.29) 



Large opening angle ^ 2). At the same time, A contributes to the multiplicity 
of the "fully open" dipole, ^ = 2. The expression p.25b|) for A(2; E) contains one 
integral: 

A(2; E) =aJ'-^ [Af{r^ ■ 2; E) - M{2- E) ] . 
Jo — V 

Expanding the difference of the multiplicity factors, 

U{'n-2-E)-U{2-E) = EdEAf{2; E) ■ Unr] + . . . , 
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we derive 



f dn TT^ 
/\{2-E)^EdEM{2-E)-as / ■ = -—a^ ■ EdEM{2- E). 

Jo ^ ^ V J-Z 

Integrating this correction over energy we finally obtain an order as correction to 
the dipole multiplicity due to emission of two soft energy ordered gluons: 

AA(2; E) = Ar°L(2; . (^1 - + . . . ) . (2.30) 

In the standard approach, this correction belongs to the coefficient function and 
originates from a "four-jet" configuration with two large-angle energy ordered soft 
gluons as additional "jets" [8]. 

2.3.3 Other observables sensitive to large-angle gluon effects 

There is a variety of observables that are sensitive to soft gluon radiation at large an- 
gles. Among them one can mention production of heavy quark pairs in jets [16], non- 
global jet observables [11] including various particle distributions and correlations 
related with particle flows in the inter-jet regions [12, 13, 17-21]. The dipole-based 
evolution equation ( p.l5| ) is well suited for taking into account single logarithmic 
correction effects due to multiple soft gluons in these observables and elsewhere. 

At the same time, it remains insufficient for building a realistic Monte Carlo event 
generator which is impossible without incorporating full parton decay probabilities, 
including hard parton splittings, and of the recoil effects to ensure energy-momentum 
conservation. 



3. Attempt to include recoil in dipole multiplication 

Generally speaking, a gluon emitted by a dipole cannot be ascribed as offspring to 
either of the two partons that form the dipole. However, in the collinear limit when 
the gluon momentum q is quasi-parallel to one of the two hard partons, say, pa, 
it is the parton a that can be said to "independently split" into two, sharing its 
momentum. Pa ^ Pa + Q- A small change in the "spectator" momentum ph 
which is necessary to compensate for the virtuality of the (a, q) pair, vanishes in the 
collinear limit. A general recoil strategy has to be formulated in such a way that 
in the collinear limit the answer reduces to the standard DGLAP parton splitting 
function Pa —>■ Pa + Q that includes hard momentum configurations, Pa ~ q. This is 
a source of a significant single- logarithmic correction, this time not from soft gluon 
radiation at large angles but from hard coUinear-enhanced emission. 

Thus, in order to properly formulate a recoil strategy, one must split the soft 
dipole radiation function Wab{q) into two pieces, wj^'^\q) and wl^l\q), each of which 



- 11 - 



incorporates the coUinear singularity when q 1 1 or q 1 1 p;,, respectively: 



{Paq){qPb) 



(3.1a) 



or, graphically, 






(3.1b) 



One can do this differently, for example by introducing [14] 

.(a), ^ _ {PaPb) uAb)r^\ - (P'^P'^) 



(3.2) 



(Paq) (Pa +Pb,q)' " (Pbq) (Pa + Pb, q) ' 

Another possibility to split the dipole is given by the "conditional probabilities" [5,6] 



w!£\q) = wM-l 



1 + 



{qPa){PbQ) - {qPb){PaQ) 
{PaPb){qQ) 



qa 



1 + 



^ab — 6 



qa 



^qb 



(3.3) 

The distribution (|3.3|) produces the exact angular ordering upon averaging over the 
azimuthal angle of the gluon momentum q around the singular direction p^: 



W. 



(a) 



ab 



■^(a-u)- 



(3.4) 



^qa 

Note that in all cases the distributions wj^l\q) (c = a,b) fall fast and essentially 
become irrelevant when the emission angle exceeds the opening angle of the parent 
dipole, > ^ab, that is away from the angular ordered kinematics. 

3.1 Recoil strategy 

Consider the elementary process 



Pa + Pb 



Pa+Pb + q- 



One ascribes definite recoil pattern separately to the two terms in (p.l|) . Thus, for 
the first term W^^^'* which is collinear singular in the direction q 1 1 p^, we choose, 
following the Catani-Seymour prescription [14] 



pi^^ = (1 - y)P,. 
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It defines the two light-like vectors Pa, Pb that represent the momenta of the parent 
partons prior to the gluon emission: 



pi") = zPa + (1 - z)yP, - kt, 
q= {I- z)Pa + zyPb + kt, 

where kt is ortogonal to Pa and P^. The on-mass-shell condition = gives 

{qPa) 



{Pa + yPb - g) = 



y 



(3.5a) 
(3.5b) 



(3.6) 



The light-cone fraction z of the parent momentum Pa carried by the final quark pa 
is given by the expression 

In the coUinear limit y ^ 0. In the soft limit both y and z 1. 

We now recast the phase space element in terms of the "parent momenta" : 

d<^ab{q) = {dpa){dpb){dq){27ry6^^\pa+Pb + q-Q) 

= {dPa){dPb){2nr5^'\Pa + Pb-Q)- (dq) X Jit\q) ^ d<^^^ ■ (dq) J^fiq). 

(3.7) 

The Jacobian of this transformation reads 

-1 



1 



{PaPb 



iPa-q,Pb) 

The gluon radiation probability we represent as a sum of two contributions: 

d^ab{q)-Wab{q) d^f^-{dq) [w!£\q) J^\q) + Wj^^ (q) J^J {q)_ 
Finally, the evolution equation in the differential form becomes 
EdEG{Pa,Pb\E) = j ^a^(qf)^(qf-Q,) 

X Y.'^SiQMbiQ)UQ)G{p'^:\q-,E)g{q,pt^-,E) - g{Pa,Pb;E) 



(3.9) 



(3.10a) 



E=uj 



In the graphic form it can be represented as 








(3.10b) 



This evolution equation is the modification of ( |2.10| ) that includes recoil as described 
above. The gluon radiation function Wab{q) = ^^'^Wab{q) is split according to 
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into the two pieces which are coUinear singular when q is parallel to pa or to pb, 
correspondingly, see [14]. The energy of the hardest real gluon qo = u = E plays the 
role of the evolution parameter. 

By introducing the Sudakov form factors this equation can be recast in an integral 
form as it has been done in the soft case in ( p.l5|) . 



g{P.,Pb;E) = SiP^,Pb;E)+ / -(^) " %f - Qo) 

J LO J Air b[Pa,Pb]uj) (3 11) 



c=a,6 



This equation generalizes (p.l5|) by taking into account the energy-momentum recoil 



and again satisfies the normalization Q{Pai Pb] E) = 1 for all u{q) = 1 (cf. (|2.11| )). 



An iterative solution of this integral equation defines a Monte Carlo process for 
generating multi-parton ensembles with account of coUinear non-enhanced single- 
logarithmic corrections due to large-angle soft gluon emission. 

Now we are going to check if the dipole recoil scheme described above is consistent 
with known analytical results concerning the coUinear resummation. The simplest 
observable of this type is an inclusive energy distribution of the final-state quark. 

3.2 The non-singlet quark energy distribution 

The non-singlet quark fragmentation function is obtained by taking the derivative 
of the generating functional Gab with respect to the quark source u{pa) and setting 
all the remaining source functions u = 1. This way we obtain the distribution in 
the momentum of the quark a accompanied by any number of gluons. If the quark 
energy is taken large, (1 — x) <^ 1, all radiated gluons are soft and the analysis 
simplifies significantly. The problem becomes essentially Abelian^ and is described 
by multiple independent radiation of soft gluons by the quark. 

3.2.1 CoUinear approximation 

We shall restrict ourselves to configurations in which all radiated gluons have small 
emission angles with respect to the quark direction. This — quasi-coUinear — ap- 
proximation is sufficient for the analysis of the anomalous dimension which accumu- 
lates coUinear ( "mass" ) singularities of the fragmentation function in all orders and 
describes the scaling violation. 

In this kinematics the measured variable x is given by 



X 



l-'^Vi, UJi = yiE, pa^xPa. (3.12) 



^Non-Abclian effects due to final state cascading of gluons reduce to making the effective coupling 

2 



in the radiation probability of a primary gluon k run with fc^ 
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What matters is the total energy 1—x carried, on average, by an ensemble of radiated 
soft gluons. 

Looking for collinear singular contributions — the terms of the type Q;sln(l — 
x) ■ InQo — "we may omit the two last terms on the r.h.s. of the equation (|3.10b| ) 
that contain the factor Wj^l\q) which is non-singular when the gluon momentum q 
becomes collinear to p^- 

One is then left with contributions involving only the singular antenna piece 






(3.13) 



Since our observable is inclusive, we may use unitarity arguments to simplify the 
equation. Indeed, by the nature of the adopted recoil strategy, production of gluons 
by the "upper" dipole (6c) in ( |3.13|) does not affect the momentum of the quark a. 
Therefore, since such emissions do not affect the measurement, they are subject to 
real-virtual cancellation and can be neglected. 

One is left to consider only gluons emitted in the lower blob involving the parton 

a: 



^dE 




(3.14) 



This equation generates multiple gluon emissions off the quark. 

By construction, successive emissions are ordered in gluon energies. And here 
comes the crucial observation: as we have discussed above, the emission angle of 
the gluon c is essentially limited from above by the aperture of the parent dipole, 
^ab, see (p. 41). Therefore, soft gluons generated by the evolution equation ( p.l4| ) 
turn out to be ordered, simultaneously, in energies and angles with respect to the 
radiating quark. Instead we know that the DGLA equation is obtaine from ordering 
only in collinear variables, disregarding the relative energies of emitted partons. The 
corresponding equation in the present dipole-formulation looks as follows. 

3.2.2 Evolution equation with dipole recoil strategy 

Since the radiated gluons are soft, we use the distribution wjj^\q) defined in (|3.3| ) 
which describes the soft part of the gluon radiation probability (splitting function). 
We may also approximately set Jal\q) = 1. 

In going from the generating functional Q{Pa, Pb', E) in ( 3.101 ) to the single par- 
ticle distribution, we express the latter as a function of the angle between the par- 
tons S,ab^ the energy scale E Ea and the energy fraction of the triggered quark: 
D{U,E,x). 
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Performing integration over the azimuthal angle (paq 

y 



D{^-E-x) = 6{l-x) + 



D 



2e 

^',yE, 



X 



y 



see (|3.4|) ) one deduces from 



(3.15a) 



D{^\yE,x) 



where ^ = ^ab is the angular aperture of the original dipole {ah) while ^' = ^ac is the 
aperture of the secondary (ac)-dipole, and for the sake of simplicity we have restricted 
ourselves to the soft part of the gluon radiation probability, dy/y. Equivalently one 
may write^ 

D{i-E-x)=5{l-x)S{E^i) 



+ 



y 



2e 



asdiyE^^'-Q, 



D{^',yE, 



X 



i-y 



(3.15b) 



In S{E^) 



as^iyE^/^' -Qo 



(3.16) 



with the Sudakov form factor given by the expression 

^'dy r^dl 

y J 2e' 

The distribution D describing real gluon emission in the integrands of ( p.l5|) is 
evaluated at ^' since a gluon produced by the dipole (ac) has an angle with respect 
to a effectively bounded from above by ^ac, see (^.4]). 

It is worth noticing that going from the equation for the generating functional 
( p.lO|) to that for the inclusive quark spectrum, we have replaced the full angle aper- 
ture ^ = C,ab by the opening angle ^' of the (ac) dipole in the virtual (subtraction) 
term as well. In so doing we followed the preceding observation according to which 
virtual radiation of gluons constituting the "wider" dipole (be) had been already can- 
celled by corresponding real emissions, without affecting the energy of the triggered 
quark. 

Equations ( |3.15| ) correctly satisfy the sum rule for the first moment: integrating 
over X one obtains 1 for the number of quarks a (to see this one uses ( 2.18| )). However, 
they are in conflict with the Poisson nature of the multiple soft gluon radiation and 
contradict the collinear resummation result leading to the DGLAP equation. 

3.2.3 Analysis of the equation 

Introducing the transverse momentum scale variable Q = Ey/E, and the correspond- 
ing integration variable q = ujy/^ we get 



D{x,Q) = 6{l-x)+ f ^ r'^as 

lo/Q y Jqo 



D 



X 



i-y 



q] - D{x,q) 



(3.17) 



^Thc y integration in ( p. 15b ) is actually limited from above by y 
D{^,Q,z) oc 9(1 -z) 



1 — X since the distribution 
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or, in terms of the standard Mellin moment representation, 

D^iQ) = 1+ f ^ r^^^a, [(1 - yf - l] (g) . (3.18) 

Jqo/q y Jqo '- 

In the region of large x, such that 1 — x ^ 1 (but 1 — x ^ Qo/Q), the essential 
Mellin moments are large, N ^ {1 — x)^^ , and cut from below the integral: 

Dr,iQ) ^ 1- r ^ r^a.D^iq). (3.19) 

Ji/N y Jqo 1 

By iterating this equation one obtains an oscillating double logarithmic series in In Q 
and In N , 

ni^)(n\ (-Q. InA^lng)" 

(Q) + .... 

Here one combinatorial factor 1/n! comes from the kinematical ordering of gluon 
energies, and the second 1/n! from the ordering of the angles with respect to the 
quark of the successively radiated gluons. This series sums up into a Bessel function 
instead of an exponent in \n.Q with the well known DGLAP "anomalous dimension" 

nWrn^ (-"s IniVlng)" (n\ n-^s\nN 

Dm [Q) , Dn[Q) ocQ ' 

where we have neglected for the sake of simplicity the running of the coupling. 

In general (whether the coupling runs or not), the anomalous dimension is 
defined as 

7Tv(as(g)) = ^l^l^^^W)- 

For as = const this gives 

Substituting this ansatz into ( |3.19| ) and evaluating the logarithmic derivative, we 
obtain an equation for the new anomalous dimension: 

Ma.) ^ -a. f dyy '"-^ = - exp(-,. IniV) ^^^^^ 

J 1/N In 

Perturbatively, when a^ln^ ^ 1, we have the correct one loop expression 

7^(as) ~ -a, InAT + C (a^in^AT) , (3.21) 

while in higher loops it starts to deviate significantly from the DGLAP answer. In 
general, introducing the variable w 

w 

= -hTiV' 
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the anomalous dimension is given by the solution of the equation 



a, In^ N 



w 



(3.22) 



The r.h.s. hits the maximum at certain w = wq, so that in the region of large moments 
N, starting from In^ = const, the anomalous dimension becomes complex valued 
and the distribution starts to oscillate. 

The origin of the failure of the energy ordering strategy can be understood al- 
ready at the level of two emitted gluons. 

3.2.4 Two gluon emission 

Consider the emission of two soft gluons pi,p2 off the parton q of the qq dipole. The 
antenna functions that potentially contribute in the collinear limit (recall that we 
keep all gluon angles with respect to the quark a = g to be small) are displayed here: 






(3.23) 



The first two graphs correspond to the splitting of the dipole (al): 



a2 



and iy«(2)^^9(e,i-e 



12 



while the third one, 

w['^{2)-.m,-ii2). 

is the relevant part of the large-aperture dipole (If^). Due to the local recoil prescrip- 
tion used, only the contribution (a) affects the momentum of the quark q. In the 
two remaining ones, (b) and (c), the gluon 2 borrows its energy-momentum from 
the gluon 1 and does not produce any quark recoil. Therefore, these contributions 
cancel against corresponding virtual corrections in the inclusive quark measurement. 
In conclusion, within the adopted recoil strategy, only the graph (a) should be kept, 
and we obtain the following phase space for the two-gluon emission: 



aq 



and 



UJ2 < OJi <^ E. 



(3.24) 



The first condition comes from the angular ordering in the graph ( |3.23| )(a), and the 
second condition — from the energy ordering of successive emissions. 

We know, however, that in order to obtain the DGLAP equation that properly 
resums collinear singular contributions, one needs to assemble angular (or, trans- 
verse momentum) ordered emissions (the first ordering), regardless to the order of 
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gluon energies (the second one). At the same time, the dipole logic is leading us 
to the double- ordered gluon ensemble, according to ( p.24|) . As a result, instead of a 
simple Poisson distribution of soft radiation (l/n!) we obtain something like a Bessel 
function distribution (l/(n!)^). 

What is missing here is actually the coherence of QCD radiation. As well known, 
a soft gluon UJ2, with UJ2 ^ uji <^ Ea, could be emitted at large angles (^a2 ^ C,ai) 
directly by the original parton g + pi ^ g. In the language of Feynman amplitudes, 
such radiation occurs as a coherent sum of the graphs (b) and (c). 

Another way to discuss this point is the following. The dipole radiation pattern 
of soft gluon 2 off the composite antenna (alq) can be represented as follows: 



Wali2) + Wig{2) 



Wali2) + Wigi2)-Wagi2) 



+ Wagi2). 



(3.25a) 



Beyond the large- A^^c approximation, the latter term (corresponding to emission of 
the gluon 2 by the quark-antiquark dipole) acquires a A^^-suppressed correction. 



Wal{2) + Wig{2)-Wag{2) 



+ CF-Wag{2). (3.25b) 



This is in line with the collinear factorization theorem. Indeed, the factor in square 
brackets is non- singular when p2 becomes parallel to the direction of the quark or 
antiquark momentum so that these collinear enhancements are contained by the 
quark dipole term Waq{2) and must be proportional to the "colour charge" (Casimir 
operator) Cp of the quark. 

The only collinear singularity of the first contribution on the r.h.s. is that of p2 
parallel to pi which suggest an interpretation of this term as independent emission 
of the gluon 2 by the previous generation (harder) gluon 1. The full two soft gluon 
emission probability takes the form 



dw{l,2) OC Cp^-Wag{ 



Wal{2) + Wi^{2)-Wag{2) 



+ Cl-Wag{l)Wag{2). (3.26) 



In the large- A'c limit, this is identical to the original dipole expression. At the same 
time, the representation (|3.26| ) suggests a different recoil strategy when it comes to 
incorporating the energy-momentum conservation. 

The first term in (|3.26| ) describes cascade gluon multiplication: qq ^ 1 followed 
by 1 — > 2. In the inclusive quark distribution this is subject to the real-virtual 
cancellation. Energy ordered configurations, UJ2 ^ i^i cancel leaving no trace; "hard" 
gluon splittings with UJ2 — uji modify one gluon emission by giving rise to as{kl 
for the radiation intensity of the primary gluon 1 by the quark. 
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The second term describes independent radiation of two gluons off the quark 
— the second order term of the Poisson distribution of primary gluons — which is 
necessary to correctly recover the DGLAP quark fragmentation function. 

Equation ( p.26|) provides a probabilistic representation for the two-gluon produc- 



tion in terms of a combination of independent and cascade gluon emissions. With 
account of hard parton splittings and of the natural energy recoil prescription, it 
served as the base for constructing the HERWIG parton event generator [5,15]. The 
corresponding probabilistic representation for production of three and four energy 
ordered gluons was developed in [6] and described in [8], which resulted in a formu- 
lation of the parton multiplication scheme based on the exact angular ordering of 
successive parton splittings. 

The difference between the two strategies becomes transparent. In the dipole 
scheme, successive gluons borrowing the energy from the quark were ordered in 
emission angles, 

while in (|3.26|) the angle of the softest gluon is not bounded, and the two emissions 
are totally independent, 

Wil\l)W^fi2), e2a<l, 6a <1. 

Let us see how this has happened, from the point of view of the gluon-gluon splitting. 

In the quasi-coUinear configuration, 6i I5 the radiation pattern of the softest 
gluon 2 consists of a "narrow" and a "wide-angle" dipoles, (al) and (Ig). When, 
within the dipole recoil strategy, the "wide-angle" antenna Wiq{2) is split into the 
pieces coUinear singular with respect to the directions of the first generation gluon 
(1) and of the antiquark, 

W,,{2) = Wil\2) + W^f{2), 

the contribution Wi^\2) is looked upon as radiation of the gluon 2 by the gluon 
1, whichever the kinematical configuration of the gluon 1 with respect to its parent 
quark (dipole factorization). At the same time, within the HERWIG logic (coUinear 
factorization), the gluon-gluon multiplication occurs only in a restricted angular 
aperture, 61 < 6a- Indeed, in the kinematical region 6a ^ 61 ^ 1 the two last 
terms in the combination 

Wal{2) + Wi^{2)-Wa^{2) 



on the r.h.s. of ( |3.25| ) cancel: 



lim 6,- 



Wlg{2) - Wag{2) 



O 



6a 

6^ 
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(while the first antenna term, Wai{2), is restricted in angle on its own). 

Physically, the gluon 2 radiated off the gluon 1 in the "wide" dipole (|3.23|) (c) at 
relatively large angles, ^21 > (,ia, is produced by the grandparent parton a (quark). 
Correspondingly, it must borrow energy from the quark rather than locally from 
the dipole (Ig) to which it belongs from the point of view of the colour connection 
topology. 

Reformulating this way the local dipole recoil prescription modifies the energy 
ordered evolution equation. The modification reduces to substituting the full dipole 
opening angle C, = C,ab for the running angle parameter ^' in the argument of the 
distribution D on the r.h.s. of (|3.15|) . As demonstrated in the Appendix, this way 
one recovers the approximate DGLAP equation for the quark fragmentation function. 



4. Conclusions 

Monte Carlo generation of QCD events is a quarter century old business, based on 
the structure of resummation of collinear enhanced Feynman diagram contributions. 
The probabilistic parton cascade picture was first established for one-particle inclu- 
sive quantities (DGLAP evolution equations for DIS structure functions and e~^e~ 
fragmentation functions, [1]), generalised via "jet-calculus" [2] and, with account of 
the soft gluon coherence [3-5,8], has laid the base for probabilistic description of the 
internal structure of parton jets and their ensembles. 

The corresponding coherent branching Monte Carlo schemes (HERWIG, in par- 
ticular) reproduce correctly both the leading double logarithmic effects (LL) in var- 
ious QCD observables as well as subleading single logarithmic collinear-enhanced 
terms. One may recall as examples NLL corrections to the mean particle multiplic- 
ity [22] and to the inclusive particle spectra from light parton [7, 23] and heavy quark 
initiated jets [24], to global event shapes in e~^e~ annihilation [25], as well as many 
other quantities sensitive to parton cascades. 

Collinear-non- enhanced ("large angle") soft gluon radiation provides significant 
NNLL corrections to global event characteristics (e.g., mean particle multiplicity); 
it is also responsible for inter-jet multiplicity and energy flows and determines the 
structure of various non-global observables [11-13]. Effects of multiple soft gluon 
radiation at large angles lie beyond the scope of the standard (collinear) approach 
and must be treated order by order in perturbation theory (while collinear enhanced 
contributions are resummed in all orders). 

An elegant expression [4] for the multiple soft gluon production probability (|2.3| ) 
is valid for arbitrary angles and offers a possibility of improving the parton picture. 
The structure of multi-gluon distribution (|2.3|) naturally suggests an interpretation in 
terms of a chain of colour connected dipoles. By choosing an "evolution time" variable 
this chain may be generated via a Markov process of successive dipole splittings. An 
evolution time can in principle be chosen differently. We have used energy ordering of 
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gluons for this purpose. The generating functional that we have constructed with the 
help of the corresponding evolution equation allows one to calculate specific effects 
due to multiple emission of soft gluons at large angles in the large- iVc approximation. 

In order to construct a realistic Monte Carlo generator for multi-parton ensembles 
it is imperative, however, to formulate an adequate recoil prescription which would 
ensure energy-momentum conservation at every successive step of the parton (dipole) 
multiplication. 

In the present paper we addressed the question, whether the "dipole factor- 
ization" extends beyond the no-recoil approximation. In other words, whether the 
splitting of a colour dipole into two can be treated independently of the prehistory 
of the system, that is locally in the evolution time (which is a necessary condition 
for constructing a Markov process). 

Having taken an inclusive energy distribution of the final quark as the simplest 
example of a coUinear sensitive observable, we have shown that a naive implementa- 
tion of the dipole recoil strategy results in violating the coUinear factorization. 

However, the existing Monte Carlo implementations of the dipole picture (for 
recent examples see for instance [26-28]) chose to order successive gluon emissions in 
(dipole centre of mass) transverse momenta rather than in (laboratory frame) gluon 
energies. Such a construction apparently prevents the problem under consideration 
from appearing in the leading coUinear approximation, that is in all orders in the 
QCD coupling corresponding to the one-loop level in the anomalous dimension. 

Indeed, as we have seen above, the problem with the quark recoil appears when 
a softer gluon (2) is emitted at a large angle than its predecessor (1). The inclu- 
sive quark energy distribution is sensitive to the region where the two gluons have 
comparable energies. Then, a softer one that carries a slightly smaller energy but is 
radiated at a large angle, acquires a larger transverse momentum, and therefore, by 
construction, is looked upon as "the first" by the A;_L-ordered MC scheme. Therefore, 
"the second" gluon in such a picture, that is the one with smaller transverse mo- 
mentum, is automatically having a smaller emission angle; in the leading logarithmic 
approximation the coherence is respected. 

The problem of an interplay between energy-momentum recoil and soft gluon 
coherence discussed in this paper gets postponed but should be kept in mind in view 
of attempts at constructing the next-to-leading order QCD MC generators. 
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A. Global recoil and correspondence with DGLAP 



We have no universal recipe under sleeve for extending the dipole gluon multiplication 
picture beyond the no-recoil approximation. However, for a simple example of the 
inclusive quark distribution discussed in Section ^ the situation is straightforward to 
cure. 

As we saw from the previous discussion, the origin of the failure lay in the 
fact that the softest gluon emitted at angles larger than the angle of the previous 
branchings did not contribute to the quark recoil. So, it suffices, by brute force, to 
permit the gluons of all generations to contribute to the quark recoil, irrespectively 
to the value of their emission angles. 



This amounts to replacing in the integrand of (|3.15b|) the distribution D eval- 
uated at ^ = ^ac with the distribution evaluated at the full opening angle ^ = ^ab- 
One then has the different equation, 

D{^,E;x) =6{l-x)S{E^0 

which actually corresponds to the DGLAP evolution. To see this, we iterate equation 
( [A.l|) to obtain the series 



Di^,E;x) = SiE^f0 



(A.2) 



where 



6en-ord = &{yn < ' ' ' < Vl < I) (A.3) 

is the product of with theta-functions which ensures the energy ordering of gluons. 
Now, due to the symmetry of the multiple integral ( |A.2| ) with respect to energy and 
angular variables, we can replace the energy ordering with the angular ordering: 

6en-ord =^ ©ang-ord = O (^n < " " " < ^1 < • (^-4) 

Now the gluons are no longer ordered in energies (the gluon does not need to be 
softer then Qi-i) but in the emission angle with respect to the quark. In the angular 
ordered form, a soft gluon radiated at a relatively large angle can be said to be 
emitted before any other (harder or softer) gluons that move in a collinear bunch 
around the quark. This is in accord with the colour coherence according to which 
a large-angle soft gluon "sees" only the total colour charge of a coUinear group of 
partons. 
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To derive the DGLAP equation we replace 



dui dyi dzi 



1 -'^Vi = zi- ■ ■ Zn (A.5a) 



where now 1 — Zi is the local fraction of energy taken away by the gluon qi from the 
parent parton (quark), while yi is the gluon energy fraction with respect to E, the 
energy of the primary parton P^: 

yi = {1- Zi)zi^i ■■■zi. (A. 5b) 

We finally obtain the evolution equation 



(A.6) 



with Q = Ey/S, and qt = yEy/^. The upper boundary of the 2;-integration, 1 — Qo/^'t, 
can be replaced by 1 in the logarithmic coUinear approximation, which results in 
the integral equation that correctly gives (the soft piece of) the DGLAP anomalous 
dimension. 



Finally, let us mention another obstacle that one faces when trying to construct 
the full realistic Monte Carlo scheme starting from the energy ordered gluon ensem- 
bles. 

In order to match the full DGLAP anomalous dimension, one has to include the 
"hard part" of the splitting function by adding to the distributions w'^^^ (q) the pieces 
vanishing at = 0. In the equation (|A.6|) this calls for the replacement 



"^'^ ~ dz,(^-lil + z^,)]. (A.7) 



However, while we could cast the leading soft distribution in ( |A.1| ) in terms of the 
local energy fraction, dyi/yi = dzi/{l — Zi), the finite term (1 + Zi) dzi cannot be 
expressed locally, via the single variable yi, see ( |A.5| ). 
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